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Abstract 

We study the {K~ ,p) reaction on nuclei with a 1 GeV/c momentum kaon beam, 
paying a special attention at the region of emitted protons having kinetic energy 
above 600 MeV, which was used to claim a deeply attractive kaon nucleus optical 
potential. Our model describes the nuclear reaction in the framework of a local den- 
sity approach and the calculations are performed following two different procedures: 
one is based on a many-body method using the Lindhard function and the other one 
is based on a Monte Carlo simulation. The simulation method offers flexibility to 
account for processes other than kaon quasi-elastic scattering, like K~~ absorption by 
one and two nucleons producing hyperons, and allows to consider final state interac- 
tions of the K~ , p and all other primary and secondary particles on their way out of 
the nucleus, as well as the weak decay of the produced hyperons into vrA''. We find a 
limited sensitivity of the cross section to the strength of the kaon optical potential. 
We also show serious drawback in the experimental set up from the requirement of 
having, together with the energetic proton, at least one charged particle detected in 
the decay counter surrounding the target, since we find that the shape of the original 
cross section is appreciably distorted, to the point of invalidating the claims made in 
the experimental paper on the strength of the kaon nucleus optical. 

1 Introduction 

The issue of the kaon interaction in the nucleus has attracted much attention in past 
years. Although from the study of kaonic atoms one knows that the iir~-nucleus potential 
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is attractive [T], the discussion centers on how attractive the potential is and if it can 
accommodate deeply bound kaonic atoms (kaonic nuclei), which could be observed in 
direct reactions. A sufficiently large attraction could even make possible the existence of 
kaon condensates in nuclei, which has been suggested in [2]. Stimulated by the success in 
reproducing the data of kaonic atoms, many works considered strongly attractive potentials 
of the order of 200 MeV at normal nuclear matter density [3-7] , or explored the dependence 
of a few observables to a wide range of depths from to 200 MeV P4T0]. More moderate 
attraction is found in similar works done in pTKTS] . Yet, all modern potentials based 
on underlying chiral dynamics of the Kaon-nucleon interaction [HHTH] lead to moderate 
potentials of the order of 60 MeV attraction at nuclear matter density. They also have a 
large imaginary part making the width of the deeply bound states much larger than the 
energy separation between the levels, which would rule out the experimental observation 
of peaks. The agreement with the data of kaonic atoms of this purely theoretical shallow 
potential is good [TH], and a fit to all data adding a small phenomenological potential to 
the theoretical one performed in [20] indicates that the best fit potential deviates at most 
by 20% from the theoretical one of |15j . 

The opposite extreme is represented by some highly attractive phenomenological po- 
tentials with about 600 MeV strength in the center of the nucleus |2T||22] . These potentials, 
leading to compressed nuclear matter of ten times nuclear matter density, met criticisms 
from [23] and more recently from [2^ , which were rebutted in [25] and followed by further 
argumentation in [26] and [27]. More recently the lightest K-nuclear system of KNN has 
also been the subject of strong debate [28143T] . 

Experimentally, the great excitement generated by peaks seen at KEK [32] and FIN- 
UDA [33l[3l], originally interpreted in terms of deeply bound kaons atoms, is receding, 
particularly after the work of [23], regarding the KEK experiment, and those of [55HS7] . 
regarding the FINUDA ones, found explanations of the experimental peaks based on con- 
ventional reactions that unavoidably occur in the process of kaon absorption. Also the 
thoughts of [38|, with opposite views to those of FINUDA in [31], and the reanalysis of the 
KEK proton spectrum from K~ absorption on ^He [32j, done in [39], where the original 
narrow peak appears much broader and is consistent with the signal seen on a heavier ^Li 
target in FINUDA [ID] , have helped to bring the discussion to more realistic terms. Never- 
theless, the possibility that the FINUDA peak of [33] could be a signal of a deeply bound 
kaon state is still defended [H]. This brief description just shows the intense activity and 
strong interest in this subject over the past few years. 

There are also claims (with very low statistical significance) of K~pp and K~ppn bound 
states from p annihilation in ^He at rest measured by OBELIX@CERN [32], as well as the 
recent claim of a K^pp bound state, seen from the pp — )■ K^X reaction by the DISTO 
experiment [^. These experimental claims are under investigation now since, before 
calling in new physics, it is important to make clear that these data cannot be explained 
with conventional mechanisms. 

In this work we focus on yet another experiment which led the authors to claim evidence 
for a very strong kaon-nucleons potential, with a depth of the order of 200 MeV ^1]. The 
experiment looks for fast protons emitted from the absorption of in flight kaons by nuclei. 
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Our aim is to show how this experiment was analyzed and which ingredients are missing. 
Throughout this work we shall see that the interpretation of the data requires a more 
thorough analysis, and with all things considered, we reach different conclusions than 
those of Ref. ^44J, in the sense that we do not find evidence for a a strongly attractive 
kaon-nucleus optical potential. 

One of the shortcomings of Ref. ^1] stems from employing the Green's function method 
[35] in a variant used in |16H38] to analyze the data and extract from there the kaon optical 
potential. The only mechanism considered in Ref. |13] for the emission of fast protons is 
the Kp — )• Kp process, taking into account the optical potential for the slow kaon in the 
final state. However, there are other mechanisms that contribute to generate fast protons, 
namely kaon absorption by one nucleon, K^N — i- vrS or K'N — t- ttA followed by decay 
of the S or the A into vrA^, or the absorption by pairs of nucleons, KNN — SA^ and 
KNN — )• AA^, followed also by similar hyperon decays. The contributions from these 
processes were also suggested in Ref. |3S]. In the present work, we take into account these 
additional mechanisms by means of a Monte Carlo simulation, while the processes involving 
KN scattering, which are dominant in this reaction, are considered in two different ways, 
one based on standard many body methods using Lindhard functions and another one 
based on a Monte Carlo simulation. The agreement of the two calculational methods gives 
us confidence to use the Monte Carlo simulation for the processes involving more than one 
step and/or one nucleon and two nucleon kaon absorption. 

2 The {K~,p) reaction in nuclei: many body approach 

We are dealing here with an inclusive reaction, where a kaon in flight hits a nucleus and 
a proton is emitted. In the present case we focus on fast protons, which could be emitted 
when the kaons are trapped in the nucleus or are rescattered with small energy. The 
reaction is inclusive in the sense that, apart from the proton observed, anything can happen 
to the nucleus. In fact many processes may take place. The original kaon can undergo 
quasi-elastic collisions with the nucleons, transferring them some energy. The kaon can be 
absorbed, either by one nucleon or by pairs of nucleons. The kaon can be trapped in a 
kaonic orbit, etc. Once a kaon has experienced a particular reaction, the final products 
also suffer their own interactions with the nucleus before, eventually, a fast proton gets out. 
Complicated as it may sound — and we shall deal with these complications in a following 
section — the evaluation of the inclusive cross section is however easy since one is looking 
for the most energetic protons and in the forward direction, which means that these protons 
must take the largest possible energy from the original kaon. In other words, if a kaon 
undergoes a quasi-elastic collision in which the proton does not fall in this narrow window, 
the event will be dismissed because this kaon, having lost a fraction of its energy, will not 
have a second chance of producing a fast proton again. Obviously, if the kaon is absorbed it 
disappears from the flux and must also be eliminated in the evaluation of further processes 
(the contribution from the absorption mechanisms will be calculated later). The other 
relevant observation is that if the final fast proton has a secondary collision it will also 
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loose energy and will not lie in the desired energy window. In practical terms this means 
that we can just care about the direct K~p — > K~p quasi-elastic reaction at a certain point 
in the nucleus and distort the initial K~ and final proton waves. These type of reactions 
have received much attention and we describe here the standard procedure of dealing with 
them [HOlET]. 

Diagrammatically, the process considered is depicted in Fig. 1, which shows the kaon 
inducing a ph excitation and remaining in the nucleus, representing what we call a quasi- 
elastic collision. The kaon transfers energy and momentum to a nucleon that is promoted 
from below to above the (local) Fermi sea. 




K-(K) 



Figure 1: Diagrammatic representation of the inclusive {K ,p) reaction. 



The kaon self-energy for the diagram of Fig. [T] in a Fermi sea is given by 



(27r)4 ^ -^^V ^ V ^ _ ^2 _ ^2^ _ n(g) ' 

where U{k — q) is the Lindhard function for ph excitation and T stands for the K^p — )■ K^p 
scattering matrix. The factor [g''^ — q"^ — f^lc ~ n(q'°, q)]'^ is the kaon propagator which 
includes its self-energy in the medium, Il{q^,q). Using the Cutkosky rules one easily 
obtains the imaginary part of nqo(A;) as in 



Then 



Hqe ^ 2^Im Hqe , (2) 

U{k -q)^ 2ie{k° - g°)Im U{k - q) , (3) 
D{q) = ,2 \ , ^ 2zg(g°)Im D{q) . (4) 



Im nqe(A:) = -2 / ^Im U{k - q)\T\Hm^^^^^ \ , (5) 
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where U is given by 

ITH, q)-2 f "^'^^ ^ ^ n{pN)[l - n{k + pN - q)] 

The fast protons with momentum k + p^ — q are the energetic ones that would be 
observed, hence the corresponding Pauh blocking factor, 1 — ra, is just unity here. iFiom 
Eq. ([H]) one obtains Im U as 

lmU = -2n f fj'l Ji,^,r ^ -n{PNnk' + E{pj,)-q'-/\-E{k + p^-q)), 

J [2Tiy E[pn) E{k + pn - q) 

(7) 

where we have introduced an energy gap between the energy of the holes and the energy 
of the particles [53li53]- We thus obtain for Im Ilqe the following equation 



xlm- 



J (27r)^ E[pn)J {2-nY E{k + pn - q , 

1 



The physical interpretation comes by recalling that 



qO=kO+E(^PM)-/:\-E{k+PM-q) 



2a;Kpt = Hqe , (9) 
Im Kpt = ;^Im Hqe , (10) 

V = -2Im Kpt = , (11) 

UJ 

where u is the kaon energy and Vopt is the /^"-nucleus optical potential. Our states are 
normalized to unity in a box of volume V. The flux of the incoming kaons is vk-/V and 
thus the K~ cross section with the nucleons of the Fermi sea of volume V is given by 

Replacing V by an integral / d^r over the nuclear density we obtain, upon the change of 
variables 

p = k + Pn - q, (13) 
where p is the outgoing proton variable, 

2 /• ,3 /• d^pN M 



a = -- d r 777^?^(pAr, r 



(27r)3 ^EipN) 
d^p ,^,2 M ^ 1 



iTr— — Im 



(27r)3' ' E{p) qO' -q"^ -ml~U{qO,q] 



, (14) 

(f=k+PN~P 



from where 
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Q=k+PN~P 



(15) 



The integral over f covers the size of the nucleus. As in previous works [35l|36] , we use a 
realistic nuclear density profile for ^^C, given by three-parameter Fermi distribution |55] . 
which reproduces elastic electron scattering data. We have also tried a more sophisticated 
nuclear density profile, which accounts for the finite range of the interaction via a folding 
procedure, and is preferred by the antiprotonic X-rays and radiochemical data [56]. Al- 
though the folded density is almost 10% lower in the center of nuclei and extends to larger 
distances, the final result of our simulation is practically unaffected by this density change. 
The integral over Pat is restricted to the hole (bound) nucleon states within the local Fermi 
momentum kpir) obtained from the nuclear density at point r. This is accounted for by 
the Pauh blocking factor n{pN,r). 

In free space, the cross section for kaon scattering off a proton in the lab frame with 
the proton emerging in the forward direction reads 



da 



lab 



k (27r) 



dn{p) 

where p is the momentum of the nucleon 

2pcM-EcM 



M 



1 
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M 



PCM 



2 p{k° + M) - E{p)k 



(16) 
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Equation ( fT6l) establishes a link between |Tp and the forward cross section, which can be 
implemented into Eq. ( 1T5|) to derive our final formula 
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TT 



02 



g 2 — m\ 



n(go,g; 



qO=k°+E(pM)-^-E{p) 



q=k+pM-p 



where we have added the distortion factor for the initial K~ and the final proton p (expo- 
nential factor in the equation), as well as the factor 6{q^) of Eq. (jlj). Taking into account 
that {a^^-j^) ~ 45 mb ~ (ap^) = cr for /('"-nucleon collisions with px- — 1 GeV/c and 
proton-nucleon collisions with protons having about 600 — 700 MeV kinetic energy, we have 
implemented a combined eikonal distortion factor as in |^I16] 



(19) 



d r — )■ / d'^re 



6 



where b is the impact parameter h"^ = + y"^. 

The backward differential cross section of the elementary process K~p — )■ K~p in the 
laboratory frame {da / drL)\s,h for incoming kaons of 1 GeV/c is taken to be 8.8 mb/sr, using 
the K~p elastic cross-section data of Ref. |57j. We note that the authors of Ref. pTj take a 
value of 3.6 mb/sr as an effective way to implement effects of Fermi motion, Pauli blocking, 
etc, which here we consider explicitly. 

3 Monte Carlo simulation 

The procedure outlined above is quite efficient to produce the cross section for the {K~ ,p) 
reaction, but obviously including only the quasi-elastic collisions K~p — )■ K~p. There 
might be other processes contributing to generate fast protons and, when this is the case, it 
becomes advisable to make a simulation of the reaction. This procedure has been developed 
in Ref. [5B] for the study of inclusive pionic reactions in nuclei and has also been applied 
to other processes, such as photon induced pion and proton emission in nuclei (SSHSl], 
electron induced proton emission [60], nucleon emission following hypernuclear decay [611 
|62] . nucleon emission following kaon absorption in nuclei [35], etc. 

As sources of fast protons we consider the quasi-elastic K~N scattering process, as well 
as the absorption of the kaon by one and two nucleons. The election of which reaction 
occurs at a certain point in the nucleus is done as usual. One chooses a step size 51 and 
calculates, by means of aip5l with i = qe, IN, 2N, the probabilities that any of the possible 
reactions happens. The values of the cross sections are discussed in Sect. |H The size of 61 
is small enough such that the sum of probabilities that any reaction occurs is reasonably 
smaller than unity. A random number from to 1 is generated and a reaction occurs if 
the number falls within the corresponding segment of length given by its probability, the 
segments being put successively in the interval [0-1]. If the random number falls outside 
the sum of all segments then this means that no reaction has been taken place and the 
kaon is allowed to proceed one further step SI. 

3.1 Quasi-elastic scattering 

We here describe how the Monte Carlo simulation treats the quasi-elastic reaction dis- 
cussed in Sect. El The general strategy is to let the kaon propagate through the nucleus 
determining, at each step 61, whether it can undergo a quasi-elastic collision, according to 
the probability Cqep^Z, where cXqe is the K~N — K~N elastic cross section. If there is a 
quasi-elastic collision at a certain point, then the initial K~ momentum and the nucleon 
momentum, randomly chosen within the Fermi boosted to their CM frame. The 

direction of the scattered momenta is determined according to the experimental cross sec- 
tion. A boost to the lab frame determines the final kaon and nucleon momenta. The event 
is kept as long as the size of the nucleon momentum is larger than the local value of kp- 
Since we take into account secondary collisions we consider the reactions K~p — )■ K~p, 
K~p — )■ K^n and K~n — )■ K~n with their corresponding cross sections. 
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Once primary nucleons are produced they are also followed through the nucleus tak- 
ing into account the probability that they collide with other nucleons, losing energy and 
changing their direction. We follow the procedure detailed in [58ll59] . 

We also follow the rescattered kaon on its way through the nucleus. In the subse- 
quent interaction process we let the kaon experience whichever reaction of the three that 
we consider (quasi-elastic, one body absorption, two body absorption) according to their 
probabilities. If the kaon remains after the collision, this procedure continues until it finally 
emerges out of the nucleus or it is absorbed by one or two nucleons. 

3.2 One body kaon absorption 

We consider the reactions K~N — )■ vrA and K~N — )• vrE, with all the possible charge com- 
binations. Once again the probability of this occurring is weighed by their corresponding 
cross sections and the directions of the vr and of the hyperons are also determined in the 
CM frame. The system is then boosted back to the lab frame, where we let the A or the S 
propagate through the nucleus, undergoing quasi-elastic collisions with the nucleons. Once 
they leave the nucleus they are allowed to decay weakly into t^N providing in this way a 
source of protons which is not negligible, as we will see. 

3.3 Two body absorption 

We also take into account the following processes: K~NN — )■ AA^ or K~NN — )■ TjN 
with all possible charge combinations. The probability per unit length for two nucleon 
absorption, iir-nn^ together with the distribution into the different possible channels, are 
discussed in Sect. HI In these reactions an energetic nucleon is produced, as well as a A 
or a S hyperon. Both the nucleon and the hyperon are followed through the nucleus as 
discussed above. Once out of the nucleus, the hyperon is let to decay weakly into nN pairs. 
Therefore, the two body absorption process provides a double source of fast protons, those 
directly produced in the absorption reactions and those coming from hyperon decays. 

3.4 Consideration of the K optical potential 

We also take into account a kaon optical potential Vopt = Re V^pt + i ImVopt, which will 
influence the kaon propagation through the nucleus, especially after a high momentum 
transfer quasi-elastic collision when the kaon will acquire a relatively low momentum. 

As discussed in the introduction, one can find in the literature quite different values 
for the real part of the potential. In the present study we vary the strength of the poten- 
tial, ReKpt, to study the sensitivity of the results: starting from — 60p/po [T^HTS] . going 
through -200p/po MeV [3HI0], and down to -600p/po MeV [2ll|22]. For the imaginary 
part of the optical potential we take ImKpt ^ — 60p/po MeV, as in the experimental 
paper [Hj and the theoretical study of [T3] . 
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In the presence of an optical potential, the kaon spectral function has the form: 



S (M ) = 1 -2M;,ImKpt .20) 

vr(M2._M|-2MxReV;pt)2 + (2MxImKpt)2' 

In the Monte Carlo simulation we implement this distribution by generating a random 
kaon mass around a central value, Mx + Re VJ^pt, which is the bare kaon mass shifted 
by the real part of the optical potential. The generated random masses lie within a certain 
extension determined by the width of the distribution Tk = — 2Im V"opt, the size of which is 
controlled by the imaginary part of the optical potential. The probability assigned to each 
value of Mk follows the Breit-Wigner distribution given by the kaon spectral function. 



3.5 Final observable 

After all the processes implemented in the Monte Carlo simulation, some particles leave 
the nucleus, and we select the events that contain a fast proton in the region of interest. 
To adapt the calculations to the experiment of [H] we keep the protons that emerge within 
an angle of 4.1 degrees in the nuclear rest frame (lab frame). For quasi-elastic scattering 
processes this would correspond to events in which the kaons emerge backwards in the KN 
CM frame and the protons are most energetic, having of the order of 500 — 700 MeV of 
kinetic energy in the lab frame. 

To facilitate comparison with experiment, the missing invariant mass of the ^^C(i^'~,p) 
reaction is converted into a binding energy of the kaon, Eb, should the process correspond 
to the trapping of a kaon in a bound state and emission of the fast proton, according to 

^J{Ek + Mi2c - Ep)^ - (Pp - Pj^Y = MiiB + Mk - Eb , (21) 

— * — # 

where Ep, Pp are the energy and momentum of the observed proton and E^, Pk are the 
energy and momentum of the initial kaon. 



3.6 Coincidence simulation 

It is very important to keep in mind that the measurements in the experiment of [H] were 
done in coincidence. The outgoing proton was measured by the KURAMA spectrometer 
in the forward direction, while another detector, the decay counter, was sandwiching the 
target. The published spectra was obtained with a requirement of having an outgoing 
proton in the KURAMA spectrometer and at least one charged particle in the decay 
counter [63] . 

Obviously, the real simulation of such a coincidence experiment is tremendously dif- 
ficult, practically impossible with high accuracy, because it would require tracing all the 
charged particles coming out from all possible scatterings and decays. Although we are 
studing many processes and following many particles in our Monte Carlo simulation, which 
is not the case in the Green function method used in the data analysis we can not 
simulate precisely the real coincidence effect. 
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What we can do is to eliminate the processes which, for sure, will not produce a 
coincidence, a procedure that we refer to as minimal coincidence requirement [68]. If 
the kaon in the first quasi-elastic scattering produces an energetic proton falling into the 
peaked region of the spectra, then the emerging kaon will be scattered backwards. In 
our Monte Carlo simulation we can select events were neither the proton nor the kaon 
will have any further reaction after such a scattering. In these cases, although there is a 
"good" outgoing proton, there are no charged particles emerging with the right direction 
with respect to the beam axis to hit a decay counter, since the K~ escapes undetected 
through the backward direction. Therefore, this type of events must be eliminated for 
comparison with the experimental spectra. 

As we will see in the next section, the main source of the energetic protons is K~p 
quasi-elastic scattering and, therefore, the minimal coincidence requirement removes a 
substantial part of the potentially "good" events changing the form of the final spectrum. 
Furthermore, events with one or two nucleon absorption or/and with several quasi-elastic 
rescatterings have a good chance of producing a charged particle that goes through the 
decay counter. Thus, the final spectrum obtained from our Monte Carlo simulations with 
minimal coincidence requirement will probably overshoot the experimental spectrum by 
an amount which will depend on the capability of the events having the given energy Eb 
of producing, apart from the corresponding energetic proton, additional charged particles 
hitting the decay counters. 

4 Input cross sections 
4.1 KN cross sections 

The elastic and inelastic two body KN cross sections for kaons of about 1 GeV/c are taken 
from the Particle Data Group (PDG) [64j. The values are the following: 

(TK-p^K-p = 21.22 mb , (TK-p^K^n = 7-15 mb , aK-n^K-n = 18.5 mb 
(Tk-p^ttOa = 4.32 mb , (Tk-p^jt+j:- = 1.76 mb 
<JK-p^-K-T.+ = 1-4 mb , cr^-p_>7rOEO = 1-58 mb 

(JK-n^TT-A = 6.35 mb , = 0.97 mb , OK-n^-K^i:- = 1-15 mb 

From the PDG we also know the total cross-sections: 

a^t.p = 51.7 mb , = 38 mb (a^-jv) = 45 mb [see Eq. (^] . 

Since these are larger than the sum of the partial cross sections that we are using explicitly, 
we define: 

aK-p^x = 14.27 mb , (TK-n~>x = 10.0 mb , 

which take care about all possible reaction channels, like K~p — > rjA, K~p — t- r/S and 
others, where no fast nucleons come out. Thus, we introduce an extra segment of length 
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(JK-p,n^x in the interval [0-1] for the Monte Carlo decision of these reactions to occur. 
If this is the case, the K~ simply disappears, since the particles produced in these reactions 
can not contribute to our observable. 



2+ 



Figure 2: The two nucleon K absorption diagrams from a kaon-exchange picture. 



4.2 Two nucleon absorption cross-sections 

The probability per unit length for two nucleon absorption is proportional to the square 
of the nucleon density: 

PK-Nn{p) = C*absP^ • 

We assume a total two body absorption rate of 20% that of one body absorption at about 
nuclear matter density, something that one can infer from data of K~ absorption in ^He [65] . 
In practice, this is implemented in the following way: 

(Pk-nn) = C^Up") = ^-^Pk-n) = 0.2«*-jv)(p) , 

where cr^x-N accounts for the total one nucleon absorption cross section and, in symmetric 
nuclear matter, it is given by: 

= + <i„ - aK-p^K-p - (TK-n^K-n)/'2 = 21.45 mb . 

Taking (p) = po/2, where po = 0.17 fm~^ is normal nuclear matter density, we obtain 

Cabs ~ 6 fm^ 
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The different partial processes that can take place in a two nucleon absorption reaction 

are: 

K~pp — )■ pA, pTP ^ nS^ 
K'pn — 7- nA, nYp , pT,~ 
K~nn — i- nS" . 



V^pj=(-60,-60) p/pg - many body approach 




Figure 3: The results of the direct many body evaluation (full line), and of the Monte Carlo 
simulation considering only the quasi-elastic scattering processes (dashed-dotted line), for 
Kpt = (-60, -60)p/po MeV. 

Ideally, their corresponding branching ratios should be obtained from relevant micro- 
scopic mechanisms, such as the kaon-exchange processes depicted in Fig. |2l There might 
be, however, other processes such as, for instance, those involving pion exchange. In the 
present exploratory work, we will consider a much simpler approach consisting of assigning 
equal probability to each of the above reactions. Noting that the chance of the kaon to 
find a pn pair is twice as large as that for pp or nn pairs, we finally assign a probability of 
3/10 for having a pS pair in the final state of K~NN absorption, 4/10 for nS, 1/10 for 
pK and 2/10 for riA. 

4.3 Nucleon and Hyperon cross-sections 

Apart from following the kaons, our calculations also need to consider the quasi-elastic 
scattering of nucleons, A and S hyperons on their way through the residual nucleus. The 
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V^pj={-200,-60) p/pg - many body approach 

V^p =(-200,-60) p/pg - Monte Carlo simuL, only QE 



250 - 




Figure 4: The same as Fig. El but for Kpt = (-200, -60)p/po MeV. 

nucleon-nucleon cross sections a^N for different momenta are taken from the parameter- 
ization of Ref. [66], as apphed also in other simulations |35|l37]. Given the uncertainties 
in the hyperon-nucleon cross sections, we may use the relation ayN = 2(T7VAr/3 based on a 
simple non-strange quark counting rule. This approximation is used for EiV scatterings. 
However, other more refined fits to experimental data also exist and, in the case of AN 
scattering, we use the parameterization of Ref. [67] , as was also done in Ref. [36] . 

5 Results and discussion 

In the first place we would like to compare the results obtained with the diagrammatic 
many-body method with those of the Monte Carlo simulation when only quasi-elastic 
scattering is considered. This is shown in Figs. [3|, H] for different optical potentials: Vi^pt = 
(— 60, — 60)p/po MeV and Kpt = (— 200, — 60)p/po MeV correspondingly. As we can see, 
the two calculations are practically identical in the region of interest. The Monte Carlo 
simulation produces slightly larger cross sections because it also takes into account the 
multiple quasi-elastic scattering processes. It is clear that these additional events, which 
contain a "good" final proton and more than one quasi-elastic collision, are rather rare. 
However, as we will see, kaon absorption mechanisms produce a substantial amount of 
energetic nucleons which need to be taken into account. 

Before showing the contributions of the new processes, let us first explore the sensitivity 
of the spectrum to the kaon optical potential. In Figs. El we show the results obtained 
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Figure 5: Results obtained using the many body method for kaon potential depths of 60 
MeV, 200 MeV, 400 MeV and 600 MeV at normal nuclear density. Experimental data are 
shown with black bars. 

with the many body method employing potential depths of 60 MeV, 200 MeV, 400 MeV and 
600 MeV at normal nuclear density. In Fig. E] we see the absolute distributions plotted 
together with the experimental spectrum [44]. Increasing the depth of the kaon optical 
potential produces an enhancement of the cross section in the bound region of kaons, as 
one might expect intuitively. The height of the theoretical distributions is much larger than 
that of the experimental cross section. We have tested that our theoretical normalization is 
correct. Indeed, if we remove the distortion of the incoming kaons and ongoing protons in 
the many body method, we obtain that the strength of the integrated cross section for the 
reaction on ^^C is six times the one of the elementary reaction, K^p K~p, at backward 
angles. The distortion implemented here reduce the nuclear cross section by about a factor 
3.5, which is also the same distortion effect obtained in ^Q\ . 

The different size of the theoretical distribution as compared to the experimental data 
is in fact showing the removal of events implemented by the concidence test applied in 
Ref. 04], demanding that some extra charged particle is detected in a decay counter sur- 
rounding the target together with the forward fast proton. It is, however, claimed in [H] 
that the required coincidence does not change the shape of the spectrum. Assuming this, 
we can rescale our calculations to give them the size of the experimental distribution, as 
is illustrated in Fig. [Hi The region of main interest corrresponds to a deep binding energy 
for the kaon (i.e. high momenta for the proton) of about 50 MeV or more to the left of the 
peak, since the quasi-elastic approach ignores many processes populating the spectrum at 
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Figure 6: The same results as in Fig. [5l but rescaled to the heigth of the experimental 
spectrum [H]. Experimental data are shown with black bars. 

low proton momenta. We observe that, to obtain a good description of the spectrum from 
—Eb ~ down to —Eb of around —100 MeV one would need optical potential depths as 
large as 400 MeV, or even 600 MeV, at normal nuclear density, not the 190 MeV claimed 
in Ref. [H], and even then we observe that the shape is not well reproduced for any of 
the potentials. 

In our opinion, these results indicate the existence of other contributions from processes 
that are not yet taken into account and/or that the assumption of an energy independent 
reduction factor due to the coincidence requirement might not be correct. These effects 
can be investigated within the Monte Carlo simulation developed in this work. In Fig. [7] 
we show the results of the Monte Carlo simulation obtained with an optical potential 
Kpt = (—60, — 60)p/po MeV, taking into account only quasi-elastic processes (dash-dotted 
line) and considering as well one nucleon and two nucleon absorption processes (solid line). 
We can see that there is a non-negligible amount of strength gained in the region of "bound 
kaons" due to the new mechanisms. Although not shown separately in the figure, we have 
seen that one nucleon absorption and multi-scatterings contribute to the region —Eb > —50 
MeV. To some extent, this strength can be simulated by the parametric background used 
in |31] . However, this is not true anymore for the two nucleon absorption processes, which 
contribute to all values of —Eb, starting from almost as low as —300 MeV. 

It is very important to keep in mind that in the spectrum of ^1] the outgoing forward 
protons were measured in coincidence with at least one charged particle in the decay coun- 
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-Eg [MeV] 

Figure 7: Calculated proton spectra with V^pt = (—60, — 60)p/po MeV, taking into account 
only quasi-elastic processes (dash-dotted line), and including all processes (solid line). 

ters surrounding the target. While a detailed simulation of these experimental conditions 
is prohibitive we can at least see their consequences by applying the minimal coincidence 
requirement. As described in Sect. 13.61 we eliminate the events that, for sure, will not 
produce a coincidence, i.e. those in which, after a primary quasi-elastic collision producing 
a fast forward proton and a backward kaon, neither particle suffer any further reaction. 
While it is clear from Fig. [7] that the main source of energetic protons in the ^^C(fC~,p) 
spectrum is the K~p quasi-elastic scattering process, many of these potentially "good" 
events will be eliminated by the the minimal coincidence requirement. As a result, the 
shape of the spectrum will change substantially, as clearly illustrated in Fig. [8] upon com- 
paring the bare spectrum obtained with a kaon potential depth of 60 MeV (solid line) with 
that obtained after the minimal coincidence cut (dashed line). The figure also shows the 
spectra corresponding to a potential depth of 200 MeV, before (dot-dashed line) and after 
the coincidence cut (dotted line). We clearly see that the sensitivity of the spectra in the 
bound region to the optical potential employed is practically lost when the coincidence 
requirement is applied. These results demonstrate the limited capability of the {K~,p) 
reaction with in-flight kaons to infer the depth of the kaon optical potential. Actually, the 
bare spectrum would be a more appropriate observable for this task. 

To further understand the effects of the coincidence requirement we introduce addi- 
tional constant suppression factors to the calculated spectrum |68] , as seen in Figs. [9] and 
[TOl Figure [9] shows our results using a shallow kaon nucleus optical potential, Kpt = 
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-Eg [MeV] 

Figure 8: Calculated ^^C(-ft'~,p) spectra for V^pt = (— 60, — 60)p/po MeV and Vopt = 
(—200, — 60)p/po MeV taking into account all contributing processes (solid and dot-dashed 
lines), and imposing the minimal coincidence requirement (dashed and dotted lines). 

(—60, — 60)p/po MeV, as obtained in chiral models. Comparing to experimental data, we 
can conclude that our results would need a reduction factor of about ~ 0.7, more or less ho- 
mogeneous in the "bound" region, —Eb < MeV, while the suppression should be weaker 
in the continuum, and basically negligible for —Eb > 50 MeV. This picture is natural from 
the physical point of view, because the spectrum to the right of the peak is populated with 
lower momentum protons. These are mostly produced in many particle final states, which 
have a better chance to score the coincidence detectors. 

However, if we look at Fig. [TOl where the calculations with a deep kaon nucleus optical 
potential of Vopt = (—200, — 60)p/po MeV are shown, we can conclude that it is much more 
difficult to obtain an overall description of the data with such a potential, even admitting 
a strong supression in the bound region and a negligible one in the continuum. 

In spite of the above described behavior, one cannot conclude that the experimental 
spectrum supports especially one potential depth over the other. However, we want to 
make clear that, in trying to reproduce the actual data, one necessarily introduces large 
uncertanties due to the experimental set up. Contrary to what it is assumed in Ref. [44], we 
have clearly seen in Fig. [8], that the spectrum shape is affected by the required coincidence. 
In fact, the distorsion of the experimental spectrum due to the coincidence requirement can 
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Figure 9: The ^"^C^K^ ,p) spectrum obtained with Vopt = (— 60, — 60)p/po MeV and the 
minimal coincidence requirement, for several reduction factors. Experimental points are 
taken from |44] . 



easily be much bigger than the difference between different potential depths, as seen by the 
sensitivity of the spectrum to the optical potential displayed in Figs. |5] and [HI Thus, the 
experiment of Ref. [H] is not appropriate for extracting information on the kaon optical 
potential. The theoretical analysis of [H] was based on the assumption that the shape of 
the spectrum does not change with the coincidence requirement. Since we have shown this 
not to be case, the conclusions obtained there do not hold. Certainly, the experimental 
data without the coincidence requirement of [H] would be a much more useful observable. 
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Coincidence simulation: U=-200 p/p^ x 1 .0 

Coincidence simulation: U=-200 p/p^ x 0.8 




-Eg [MeV] 

Figure 10: The same as Fig. El but for Kpt = (-200, -60)p/po MeV. 
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